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ABSTRACT 



Analyses and interpretation of surface layer and synoptic-scale data 
obtained in the Northeast Atlantic were performed to obtain descriptions 
of the evaporation duct and associated atmospheric and oceanic synoptic 
features. The surface layer data were quite unique because they were ob- 
tained from high quality measurements from ships spatially separated in a 
fixed array. Magnitudes and horizontal homogeneity of duct heights were 
compared to air-mass trajectories and weather patterns. The mean duct 
height was 4- 5m and higher values (8- 13m) occurred with North and North- 
west trajectories. During times when the heights were 8-13m, horizontal 
homogeneity occurred 35% of the time with a maximum duration of 12 hours. 
This was established on the basis of a two hundred kilometer separation 
between locations of duct height estimates. 

Significant features of the duct height and its temporal and spatial 
variations were related to synoptic scale descriptions. Satellite imagery 
used in conjunction with point observations appears to provide the most 
useful information in describing the intensity and areal distribution of 
the evaporation duct. 
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I . INTRODUCTION 



Over the surface of the ocean there exists a phenomenon described 
as the evaporation duct. The military significance of its existence 
has been known for many years. Attempts have been made to quantify 
the evaporation duct height and determine its effect on electro- 
magnetic propagation. Assessment of evaporation duct characteris.tics 
in tactical military applications should only require routine ship- 
board measurements and simple calculations; hence, a great deal of 
study has been performed to parameterize the duct on the basis of 
simple observational measurements of surface and near-surface para- 
meters. Also, as the use of low-flying supersonic aircraft and 
missiles become increasing threats to Naval surface combatants , timely 
and precise descriptions of the nature of the evaporation duct will be 
required to permit its exploitation. 

It has been expected that in an area in the high latitudes, well 
away from large land-mass influences, and where frequent migratory 
weather systems disturb the nature of the atmosphere near the surface, 
the height of the evaporation duct is likely to be small and its 
occurrence highly localized. Past research has been primarily limited 
to regions where ducts are known to be both strong (high) and per- 
sistent (e.g. tropics, southeast portion of subtropical high, coastal 
regions of warm offshore winds) with relatively little examination of 
higher latitude, open oceanic regions. This study will examine the 
occurrence and nature of the evaporation duct over the North Atlantic 
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Ocean subject to highly varying atmospheric conditions. This 
examination will focus on determining significant ducts and their 
1) frequency of occurrence and duration, 2) temporal and spatial 
nature, and 3) sensitivity to atmospheric conditions. 

The results reported here were made utilizing a unique set of data 
acquired from the Joint Air-Sea Interaction Project (JASIN) conducted 
in the Northeast .Atlantic during the late summer of 1978. The results 
were analyzed and comparisons were made bet\\?een several ships, syn- 
optic analyses and high resolution satellite imagery. The resultant 
examination yielded a description of an evaporation duct character- 
istic of the region as it was affected by frequently changing weather 
patterns . 

This study was on data collected from four research vessels 
engaged in the JASIN experiment from mid- July to early September, 1978 
under open ocean conditions. The quantity of data and the large area 
of study permitted an examination to determine the nature of the duct 
and the atmospheric disturbances most likely to influence it. An 
important aspect of the experiment was the positioning of four ships in 
a triangle (180 x 180 x 200 km) during periods of intensive meteorologi- 
cal measurements. This constant spatial separation permitted accurate 
assessments of spatial as well as temporal variations over a large 
horizontal region as weather systems moved through the area. The data 
were used with expressions presented by Fairall et al [19781 to compute 
the associated oceanic evaporation duct height. 
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II. BACKGROUND 



It has been long recognized that distributions of temperature and 
humidity in the lower atmosphere are responsible for anomalous electro- 
magnetic (EM) wave propagation ( Katzin, 1947 and Kerr, 1951 ] result- 
ing in either enhancement or degradation of Naval sensors and weapons 
systems. In particular, propagation between terminals reasonably close 
to the sea surface has been observed to exhibit significant variation 
of signal strength at ranges far in excess of the radar horizon. 

(Fig. 1) Such anomalous propagation effects were first seriously 
studied during and immediately following World War II. In recent years, 
advances in mixed- layer theory have led to numerous computational 
approaches for real-time assessment of atmospheric refractive effects. 

Often relatively cool moist marine air extends vertically from 
the ocean surface to altitudes of a few hundred meters. Above these 
heights the air becomes much warmer and drier for a variety of reasons. 
Therefore, a transition layer exists in which the air ten^erature in- 
creases and moisture rapidly decreases with increasing height. The 
resulting gradients can cause the EM refractive index to decrease 
rapidly enough to refract low incident angle radio waves downward 
relative to the earth and create ducting conditions. If the decrease 
in refractive index is great enough a surface-based duct occurs which 
will trap surface-originated radar signals and give greatly extended 
ranges for detection of near surface targets. 
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Figure la. Radar wave path under "standard” atmospheric conditions. 

Note path curves downward but at a rate less than the earth’s curvature. 
Beyond- the-horizon target detection is not possible. 




Figure lb. Radar wave path under ducting conditions. Path curves down- 
ward at a rate exceeding the earth's curvature resulting in beyond-the- 
horizon target detection. 
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The air adjacent to the ocean is saturated and the relative 
humidity is near 100%, decreasing rapidly in the first few meters to 
ambient values which depend on varying meteorological conditions. The 
decrease of humidity causes the modified refractive index (M) to 
decrease with height initially. At greater heights, the humidity 
gradient decreases which will cause M to reach a minimum and thereafter 
increase with height (Fig. 2). The height at which M reaches a 
minimum value is called the evaporation duct height (Z*) and is a 
measure of the strength of the evaporation duct. 





Figure 2. Relative humidity and modified refractivity M versus 
altitude for an evaporation duct. 
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The probability of occurrence of duct heights large enough to 
cause beyond- the-horizon detection capability for a particular radar 
varies according to geographic location, season and the time of day. 
The strongest ducts are generally found at tropical latitudes, during 
the summer season, and during daylight hours. 
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III. EVAPORATION DUCT MODEL 



Flux- profile relationships in the atmospheric surface layer 
(Z<50m) above the ocean surface have been studied extensively and are 
well understood. Bulk methods can be used to determine surface flux 
scaling parameters from routine shipboard measurements. These scaling 
parameters can be applied to refractive equations to determine 1 * . 

The method requires only mean meteorological conditions as input para- 
meters and is vvell suited to assessments using routine shipboard obser- 
vations . 

A. EM PROPAGATION CONSIDERATIONS 

The refractivity N (in N units) of the atmosphere for EM waves is 
given by Battan [1959] 

N = 77.6 P/T + 3.73x10^ e/T^ (1) 

where P is the atmospheric pressure, T is the absolute temperature 
and e is the partial pressure of xvater vapor. An expression for the 
vertical gradient in N is given by 

dN/dZ = Cl 6P/6Z + C2 5q/5Z + C3 6T/6Z (2) 

where q is the water vapor mixing ratio (q=.625e/P in gm/kg, P in mb, 

T in °K, and Cl=.3, C2=7.2 and C3=-1.3). A negative gradient in refrac- 
tivity at height Z causes a horizontally propagating EM wave to be bent 
to^vard the surface and if the gradient is large enough, the amount of 
bending exceeds the curvature of the earth and the wave becomes ducted. 
Since magnitudes of specific humidity gradients decrease with increas- 
ing height in the near-surface layer, the gradient is no longer strong 
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enough for ducting at some height I*. This height is termed the evapo- 
rative duct height aind is defined by the critical gradient necessary 
for EM trapping. 

dN/dZ = -.157 (m‘^) 

using Equation (2) and noting dP/dZ= -pq (the hydrostatic balance for 
the atmosphere) , 1* is then defined by the following equation. 

-.157 = -.032 + C2 6q/6Z + C3 6T/6Z (3) 

Since 6q/6Z and 6T/6Z are height dependent (approximately 1/Z) there 
exists a value of 1=1* such that Equation (3) is satisfied. The bulk 
method allows one to relate 6q/6Z and 6T/6Z, which are needed to deter- 
mine 1* , to the surface layer fluxes and, hence, to mean parameters. 

B. BULK .METHOD 

The theoretical framework for determining the evaporation duct 
height comes from knowledge of surface layer scaling parameters using 
the ivfonin-Obukhov similarity (MOS) theory. The theory relates the sur- 
face layer profiles of temperature, wind speed, water vapor and turbu- 
lence to the surface fluxes of momentum, sensible and latent heat. 
Vertical gradients of wind speed (U) , potential tenperature (T) or water 
vapor mixing ratio (q) are given by the following general expression 
6X/6Z = (XVukZ) p^(Z/L) 

IVhere X = U, T or q, k = Von Karman's constant (0.35), a is a dif- 
fusivity constant, X* is the scaling parameter and p(Z/L) is the >05 
stability correction. The surface fluxes define the scaling parameters 
within the flux-profile relationships. Scaling parameters for 
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velocity, potential tenperature and specific humidity are related to 
the normalized fluxes by: 

U* = (-u'w' )'^ (m/s) 

T* = (-T^')/U* (°K) 

q* = (-^')/U* (g/Kg) 

where u'w', T'w' and q'w' are the mean turbulent fluxes of momentum, heat 
and moisture respectively. A scaling length, the iVfonin- Obukhov length 
(L) , is used to account for stability influences. Stability magnitudes 
are based on the height scaled by L, i.e. Z/L. Negative and positive 
Z/L values corresponded to unstable and stable conditions. L is defined 
by the scaling parameters as 

L = TU*^ /(Kg(T* + .61q*T)) (5) 

The bulk method is based on integral forms of Equation (4) . Integrating 
Equation (4) from Z = Z^^ to some reference height, Z, an expression 
for the value of a parameter at the measurement height is obtained 
[Businger, 1973] • 

X(Z) = X(Z^^) + XVaKCln Z/Z^^ - i|;(Z/L)) (6) 

where }p(Z/L) is the profile stability function and subscript 
X = U, T, or q. (Note: a and ^ depend on X) 

Utilizing Equation (6) , a method of estimating surfaces fluxes in- 
volves the scaling parameters and roughness lengths, Z^^^. Appropriate 
values of tenperature, relative humidity and wind speed are obtained 
at both the sea surface and at some reference height, Z. Equation (6) 
can be used to determine scaling as a function of the differences. 
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= (kU) / (In (Z/Z^) - ij;i (Z/L)) 

T* = (T-Ts) aK/ (ln(Z/Z^^) - 1^2 (Z/L)) 
q* = (q-qs) aic/ (ln(Z/Z^^) -iPz (Z/L)) 
a was set equal to 1.35 for T and q from Businger et al [1971] . The 
quantities Z^ and Z^.p are roughness lengths for velocity and tenpera- 
ture profiles. It is convenient to rewrite these equations in the 
generalized drag coefficient form 



X* = (X(Z) - X(0)) 



(7) 



The neutral stability drag coefficient can be defined as 



= a</ln(Z/Z^^) (8) 

then the drag coefficient becomes 

/ (1 - ((C^S /a<) (ip (Z/L ))) 2 (9) 

In this form, Z^^ is related to the neutral drag coefficient by 

= cc</ ln(Z/Z^^) (10) 

where roughness length 

2 

“oX ^ ^ ^ (11) 

Equation (S) can be rewritten to give the ^DS stability parameter, 
Z/L, as 

(1 - /< (>1,1 (Z/L ))) 2 

Z/L = (Z/L)^ 



( 12 ) 



(1 - (C.^^) /a< (iPz (Z/L))) 



where 



C ^ 

KgZ ^TN 



((T(Z)-T(O)) + .18(q(Z)-q(0))) 



(13) 



(Z/L)^= — 



Sn 



U(Z) 
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Note that (Z/L)^ is an initial estimate of stability in terms of 
differences and neutral drag coefficient. Final values of Z/L and 
scaling parameters can be found from a single iteration process. 
Stability correction fmctions used were: 

Velocity profile: 

i|;i(Z/L) = 2 In ((1 + X)/2) + In ((1 + X^/2) 



-2 ARCTAN (X) + n/2 Z/L<0 

X = (1 -15 Z/L)^ 

iJ;i(Z/L) = -4. 7 (Z/L) Z/L>0 

Tenperature profile: 

ij; 2 (Z/L) = 2 In ((1 + X) /2) Z/L<0 

X = (1 -9 Hh'p- 

ijjzCZ/L) = -6. 5 (Z/L) Z/L>0 



The velocity drag coefficient, in a near- neutral atmosphere 

was calculated from the wind-speed dependent formulation given by 
Kondo [1975 ] . A typical value of at Z=10 m was 1.3x10“^ m/s 
which yielded Z^=6xl0 '* m. Wind speeds and drag coefficients are 
given in Table I. 
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TABLE I 



VERSUS WIND SPEED AT 10m [ KONDO, 1975 



U(m/s) 






.3-2.2 
2. 2-5.0 

5. 0- 8.0 

8.0- 25.0 



1 . 08 xU'-^^ 

.77 + .086 X U 
.87 + .067 X U 
1.2 + .025 X U 



The temperature and moisture drag coefficients, and at Z=10m 

■ 3 



■'TN ' — • '"qN’ 

were estimated as 1.3x10 MDavidson et al, 1978]* For ct.j.=1.35 
[Businger, et al, 1971 ], a roughness length of Z^.p=2 .0x10" ^ m, was 
obtained, which is assumed to be independent of wind speed. 



C. Z* EQUATION IN MOS FORM 

In the surface layer, the gradient of N can be written in terms of 
T* and q* as 

dN/dZ = -.032 = q(Z/L)/CaKZ) (7.2q* - 1.3T*) (14) 

which arises from Equation (3) . The duct height , Z* , is obtained from 
Equation (14) by setting dN/dZ = -.157, the critical gradient for ducting. 
-(7.2q* - 1.3T*) 

Z* = q(ZVL) (15) 

aK(.125) 

Given q*, T* and L, Equation (15) was solved iteratively to obtain Z*. 

No ducting occurs unless 
7.2q* -1.3T*< 0 
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IV. DATA ACQUISITION 



The data used in this analysis were unique in that they represented 
accurate measurements of the atmospheric surface boundary layer on a 
horizontal scale of 1-200 km. All data were gathered from scientific 
research vessels whose primary cruise objectives were to observe, dis- 
tinguish and qtiantify fluxes and mixing in atmospheric and oceanic 
boundar>’ layers. The following are descriptions of the measurement 
platforms, their locations, the data accuracy and the intercomparison 
of data between ships. 

A, PLATFOR^IS AND LOCATION 

Data for this analysis were taken during the JASIN experiment in the 
summer of 1978. Data from four ships involved in this experiment were 
selected on the basis of two objectives: 

1. Examine as much of the period as possible. 

2. Examine variation over nearly fixed spatial separation. 

The four ships selected were the METEOR, the JOHN MURRAY, the GARDLINE 
ENDURER and the HECLAo Figure 3 illustrates the profile and sensor loca- 
tions for each ship. Hourly observations were made of air temperature 
(T) , wet bulb temperature (Tw) , sea (bucket) surface temperature (Ts) 
and wind speed (U). The triangle configuration of the ship locations, 
depicted in Figrore 4, provided quite good areal coverage. Tlie ships 
were positioned approximately 200 km apart at each comer of the 
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R.V. \ETEOR 



A = 23ni 
B = 11m 




©' 




Figure 3, Ship profiles and sensor locations: (A)- U, 
(B) - T,Tw and (C) - Ts. 
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triangle and one ship was located at the center. Nominal positions 
were: 



MA - 59.00 N 12.5 W 
MB - 60.25 N 10.5 W 
MC - 60.25 N 14.5 W 
CENTER - 59.77 N 12.5 W 



B. DATA COLLECTION 

The subject observations were recorded on log sheets and plots of 
three-hourly values were prepared by JASIN participants. The plots 
from the GARDLINE ENDURER, JOHN MURRAY and HECLA, were the primary 
source of data used in this study. Air temperature, wet bulb tempera- 
ture, sea surface temperature were accurate to 0.2 °C, pressure to 
1 mb and wind speed to 1 kt. .Additionally, data used from the METEOR 
were available in computer printout. 

C. DATA QUALITY AND INTERCOMPMISON 

During the period 15-22 July, on 19 and 20 August, and on 
5 September, the four ships were positioned for an intercomparison of 
meteorological measurements. Hourly (half-hourly during formal inter- 
comparisons) values of wind speed, atmospheric pressure, dry and wet 
bulb air teiiperatures and sea surface tenperature were compared between 
vessels approximately one kilometer apart. These data were initially 
evaluated for coding errors and inconsistencies resulting from position 
changes. Mean and standard deviations of each variable were calculated 
for each participating vessel. Time- series of inter-platform sensor 
disagreements were produced from mathematical differences between 
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obser\’’ed variables . Interplatform means and standard deviations of 
these differences were calculated and, assuming the disagreements to be 
normally distributed, the 90% confidence interval of the estimates of 
differences were calculated. 

Results derived from this procedure were time records of inter- 
platform sensor differences. Constant disagreement throughout JASIN 
implied the sensor performed in a consistent manner and it remained 
only to explain these errors in terms of knonn disparities, such as 
non-uniform sensor heights or instrument bias. A variation of the 
interplatform corrections required further analyses either in terms of 
atmospheric stability or in terms of trends of sensor drifts. The 
latter were established by regression techniques. Inconsistencies 
unexplained or uncorrectod by these considerations remained as errors 
in the four-ship system. 

The I'^ETEOR data were selected as the standard to which data from the 
other vessels were adjusted. However, the METEOR'S sea surface tempera- 
ture values were corrected with data from the other three ships . The 
NETEOR was chosen because it was involved in more intercoraparisbns , 
because it had a preferred method for station-keeping, and because of 
its proximity to other JASIN ships and buoys . The corrections accord- 
ing to measured parameters and individual ships are summarized below 

[Macklin and Guymer, 1980] • 

1. Wind speed (m/s) 

NETEOR = METEOR 
HECLA = 0.89 (HECLA) + 0.3 
MURRAY = MURRAY 
ENDURER = ENDURER 
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2 . Pressure (mb) 



METEOR = METEOR 

HECLA = HECLA + .08 To 0600, July 28 
= HECLA + .40 From 0900, July 28 
MURRAY = 1.00640 (MURRAY) - 6.23 
ENDURER = ENDURER + .38 To 1700, August 19 

= ENDURER + .38 After 1700, August 19 (P less than 1000) 

= ENDURER + .48 After 1700, August 19 (p greater than 1000) 

3. Dry Bulb Temperature (®C) 

NETEOR = METEOR 
HECLA = HECIA - 0.3 
MURRAY = MURRAY - 0.5 
ENDURER = ENDURER - 0.5 

4. Wet Bulb Temperature (°C) 

NETEOR = NETEOR 
HECLA = HECLA - 0.3 
MURRAY = MURRAY - 0.3 
ENDURER = ENDURER - 0.1 

5. Sea Bucket Temperature (°C) 

METEOR = NETEOR + 0.3 To 1300, July 27 

= METEOR 1300, July 27 to 0100, August 3 

= METEOR - 0.4 0100, August 3 to 1200, August 9 
= METEOR From 1200, August 9 

HECU = HECLA 
MURRAY = MURRAY 
ENDURER = ENEURER 
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V. ANALYSIS PROCEDURES 



Data analysis procedures were designed to determine 1) the repre- 
sentative Z* values for the summer regime, 2) meteorological factors 
which most affected Z*, 3) the frequency of occurrence and duration of 
ducts, and 4) the extent of horizontal homogeneity of the duct. To 
achieve this, the data were first corrected, then used as input to 
bulk formulae to compute Z*. 

A data file was constructed of the most representative Z* values 
and corrected Ts-T, U and RH at three-hour intervals. Winds were 
reduced from anemometer height to temperature sensor heights using the 
stability- corrected logarithmic wind profile [ Businger et al, 1971 ]. 
The periods examined and sensor heights are summarized in Table II. 
(Dates are in year-month-day-hour format) 

TABLE II 



SHIP DATA PERIODS AND MEASUREMENT HEIGHTS 



SHIP 


STARTING 


ENDING 


ANE^DMETER HGT (m) 


Z(m) 


MURRAY 


78071312 


78072103 


11 


6 


METEOR 


78072106 


78080912 


23 


11 


MURRAY 


78080915 


78081012 


11 


6 


METEOR 


78081709 


78090506 


23 


11 


ENDURER 


78090509 


78090606 


15 


4 



Statistical and time-series analyses of Z* were performed along 
with subjective analyses of influencing meteorological conditions. The 
latter were performed with individual observations, synoptic weather 
maps , and satellite imagery (visual and infrared) . These descriptions 
were correlated to the time-series of Z* values. Joint probability 
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density distribution and conditional mean function results were chosen 
to examine correlations of at least two weather parameters to Z* 
values. Finally, con^jarisons were made of Z* and all atmospheric 
measurements to assess the horizontal homogeneity of the duct for 
periods when ships were spatially separated and at fixed locations. 

VI. PRESENTATION AND INTERPRETATION OF RESULTS 

A. GENERAL ANALYSIS OF SUMNER PERIOD (13 July - 6 September) 

Figure 5 is a time- series comparison of evaporation duct height 
(Z*), air-sea temperature difference (Ts-T), wind speed (U) , and 
relative humidity (RH) for the period of study, 13 July to 6 September, 
with the exception of the interval between 1200 CNff, 10 August and 0900 
GMT, 17 August. This time series was formulated from at least one ship 
located in the JASIN array and provided a near -continuous analysis of 
the evaporation duct height (Z*) at three-hour intervals . General 
conditions during the period will be described in the following sections. 
1. Synoptic Conditions in JASIN Area 

Daily weather maps for the Northeast Atlantic are shown in 
Appendix A. Julian day and date (day of the month) are shown on each 
map (see also Table A-I). Studies of the daily surface pressure charts 
and hourly observations [ Guymer and Taylor, 1978 ] showed that there 
were essentially four synoptic situations affecting the JASIN area 
during the entire period. 
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Figure 5. Time-Series (13 Jul-6 Sep) analysis of Z*, Ts-T, U and RJl with synoptic 
weather types, special atmospheric events and mean values indicated. 



a. Moist Anticyclonic (A) 

The periods of 14-18 July and 22-30 August were affected by 
moist westerlies or northwesterlies on the northern flank of an anti- 
cyclone to the west of Ireland. Weak frontal waves were embedded in 
this flow and occasional incursions of drier air took place. Winds 
varied between 3 and 12 m/s. 

b. Mobile Westerly Conditions (B) 

The periods of 19-21 July and 30 August - 2 September were 
characterized by depressions moving through the area bringing cold and 
dry air. The disturbances were very weak. 

c. Cyclonic Southwesterly (C) 

During the periods 22-29 July; 11-21 August and 2-5 Sep- 
tember the depressions became slow moving to the south of Iceland and 
fronts or troughs from the southwest affected the region. Generally, 
the fronts were occluded and slow moving by the time they reached the 
area, with only moderate winds. However, winds which were the strong- 
est experienced during the entire JASIN period reached speeds in excess 
of 15 m/s on 22/26 July, 18/20 August and 5 September. The air was 
significantly warmer than the sea. 

d. Blocked (D) 

During the period 30 July - 10 August a blocking low formed 
over northwest France and moved very slowly northeast; to the west a 
ridge persisted. Low level northeasterlies were colder than the sea 
surface but convection was inhibited by v^^arm air aloft associated with 
weak, slow moving occlusions which also produced occasional drizzle. 
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Towards the end of the period, the ridge intensified and moved eastward, 
e. Special Events 

Some interesting atmospheric occurrences were noted as 
events (refer to events marked on Figure 5) and are listed below. 

1. Relatively cold, dry air (14, 20 July; overnight 1-2 August; 

7-9, 29, 51 August) . 

2. Significant moisture flux in near-neutral conditions (30 July; 

23-25 August; 1 September). 

3. A warm frontal wave on 25 August and a cold front preceded by 
convective instability on 28 .August. Surface fluxes were very small. 

4. Winds in excess of 15 m/s (22 July; 18, 19 August). 

5. Incursions of warm, dry air (14, 18, 25, 28, 29 July). 

The following are descriptions of events which caused warm, dry air 
(not a normal occurrence over the open ocean) to occur. 

(1) Afternoon Heating . On 14 July, both the >IURRAY and 
HECLA reported increasing- air ten^erature but MURRAY reported three 
consecutive wet bulb temperature readings markedly lower than HECLA. 

On 18 July, MURRAY and METEOR had similar air temperature readings 
while HECLA was 1.4®C lower. All three measured similar values for 
wet bulb temperature. On both occasions, winds were light under 
prevailing high pressure. 

(2) Subsidence and Continental Trajectory . On 25 July, 

a weak ridge moved over the area as a low west of Ireland moved north- 
east. Winds backed from the northwest to the east -southeast. Warm, dry 
air, a result of subsidence from earlier ridging, flowed from the 
British Isles and influenced the area for six hours before the approaching 
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occlusion became dominant. 

(3) Short Wave Troughs . Throughout 28 and 29 July, a long 
wave trough was oriented north- south from Iceland. Weak short wave 
impulses on both days provided surges of dry air, each lasting nine hours, 
while generally maintaining warm southerly flow. 

2. Sea Surface Temperature (TS) 

Surface water characteristics were associated with the North 
Atlantic Drift, an extension of the warm waters of the Gulf Stream. The 
mean water temperature in the vicinity of the JASIN triangle was 12.2°C. 
The general areal distribution of the surface temperature for the 
sampled periods is depicted in Figure 6 [Liu 8 Katsaros, 1978]. 

The Northeast comer of the triangle was always the coldest position 
(11.2°C). The Northwest and South comers of the triangle were always 
the wamest (between 12 and 13°C) . A large tongue of cool water 
(11.8°C) from the north, intruded south and south- southeast through 
the center of the triangle. It persisted for the entire period result- 
ing in minimal temperature variations at the center position (about 
0.2°C). The most pronounced Ts gradients occurred along a line betiveen 
the NW and NE stations. 

3. Representative Z* Values 

A statistical analysis of observed parameters and calculated 
values from 714 independent observations was performed to deteimine a 
representative Z* value (Table III). The mean evaporation duct height 
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Figure 6, Composite Sea Surface 
Temperature Analysis 



was 4.2 meters with a standard deviation of 3.9 meters. This height 
appears low when compared with the climatology of the region. Sweet 
[1979] determined an* evaporation duct climatology at Ocean Station 
-ALFA (62N/33W), 1100 km upstream of the JASIN array, and Ocean 
Station INDIA (58N/19W) , 400 km southwest of the JASIN array. The 
following statistics are a comparison of the JASIN evaporation duct, 
using the method previously described and climatology which used bulk 
aerodynamic formulae and a method originated by Jeske [1971] and modi- 
fied by Hitney [1975 ]. 





ALFA 


JASIN 


INDIA 


MEAN (m) 


5.8 


4.2 


7 


^EDIAN (m) 


5.7 


3.8 


7 


IQR (m) 


3. 4-8. 3 


1.9-5. 7 


4.1-9 


# OBS 


7058 


714 


6758 



The interquartile range (IQR) represents the lower and upper bounds for 
the middle duct heights (251 of the ducts have heights below the lower 
number and 75% have heights below the upper number) and can be viewed 
as an indicator of the spread or variation of the duct height. 

TABLE III 

STATISTICS FOR JASIN PERIOD 13 JULY - 6 SEP 1978 
(714 OBSERVATIONS) 

TS T TS-T U RH Z/L T* U* q* 1* 
(°C) (°C) (°C) (m/s) (%) CQ (m/s) (g/kg) (m) 

MEAN 12.2 11.8 .5 6.6 90 -.1 -.02 .27 -.05 4.2 

SIGMA .6 .8 .9 2.9 8 1.0 .9 .13 .03 2.9 
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4. Atmospheric Factors Affecting 1* 

Since 1* is a conputational artifact determined by four meteor- 
ological parameters, Z* (Ts-T, Rh, U, Z/L) , multivariable analyses 
were used to determine the degree to which individual variables affected 
duct strength. Joint probability density distributions were computed 
for each of the parametric interrelationships to examine the pattern of 
atmospheric factors as they are related to the occurrence of 1* values. 

Based on probability theory [Batchelor, 1953] and utilizing a 
technique described by Holland [ 19681 , the distribution of 1* was ana- 
lyzed statistically against the distribution occurrence of a combination 
of meteorological variable such as Ts-T versus U, (Figures 7-12). 
Individual moments and distributions for each variable were also computed 
in the analysis. The following data were tabulated: 

1 . Mean 

2 . Variance 

3. Standard Deviation (Sigma) 

4. Skewness 

5 . Kurtosis 

6. Dimensional Class Limits by Sigma/4 Intervals From -4.0 Sigma 
to +4.0 Sigma 

7. Density of Observations for Each Occurrence by Sigma/4 Class 
Interval 

Individual class and joint class intervals were computed for 
each observation set. Two indices were computed for each data pair 
representing its position in an array (18 x 18) of joint class inter- 
vals (Sigma/2 x Sigma/2) centered on the mean values of the two 
variables. Class intervals were determined by dividing the standard 
value of one of the joint variables by 0.5. Pairs of independent 
variables were then used to compute joint probability density functions 
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Fi^re 7. Conditional mean (plotted numbers) of 1 * for 
joint -frequency distribution (delineated with solid labelled 
isolines) of wind speed, (U) and air- sea temperature differ- 
ence, (Ts-T) . Integers along the right and bottom edges are 
sums of occurrences in respective rows or columns. The total 
sum is 714 but may be less if the compared variables are 
truncated when out of range. 
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Figure 8 (top) . JDF of Z/L and Figure 9 (bottom) JDF 
of RH and Ts-T to CMF of 1 * . 
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Figure 10. JDF of wind speed and relative humidity 
to CMF of Z*. 
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Figure 11. JDF of relative humidity and stability 
to CMF of Z*. 
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Figure 12. JDF wind speed and stability to OIF of Z*. 
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(JDF) . Similarly, a conditional mean function (OIF) , defined as the 
average values of Z* for all observation times when the two independent 
variables had values within the joint intervals, was defined for each 
joint class interval. A subjective analysis of the JDF's and the OIF’s 
was performed to distinguish relationships as being from computational 
artifacts or from accidental occurrences. 

a. Air-Sea Temperature Difference (Ts-T) : (Figures 7, 8 and 9) 
The air temperature averaged colder than the sea sur- 

face temperature. Z*(Ts) calculations were largely influenced by the 
observing ship's location in the JASIN area (Fig. 6). Air temperature 
variations were primarily influenced by migratory weather systems. The 
most notable positive deviation of Z* occurred when Ts-T was less than 
zero (-.5°C), with low relative humidity (RH<85%) and stable atmospher- 
ic conditions (Z/i> 0) . However, under these conditions duct strength 
was extremely sensitive to humidity and stability changes. If RH in- 
creased Sigma (RH) or Z/L decreased Sigma (Z/L) then the evaporation duct 
virtually vanished. Generally, Ts-T had little effect on Z* when 
RH>90%. As Ts-T approached zero or became positive, more positive Z* 
deviations occurred for a wider range of RH e.g. (50-90%). The poorest 
ducting conditions occurred when T exceeded Ts more than 0.5®C and RH 
was in the high 90% ’s. In near-neutral conditions, as Ts-T increased 
Z* tended to increase fractionally. Generally then, Ts-T was more of a 
determinant than stability, except in the special situation of warm, 
moist air over cooler water. 
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b. Relative Humidity: (Figures 9, 10 and 11) 

Average humidity was 90% and ranged from 49% to 100%. 
Significant changes in humidity were directly due to the passage of 
large-scale weather systems. As relative humidity dropped to 80%, 
marked increases in Z* occurred. The largest Z* values occurred when 
dry (RH<80%) and warm air overlay cool water (Ts-T <-.5'’C). For the 
relative humidity range (50-90%), larger than average I* values 
occurred when the air was 1.0°C cooler than the sea surface. Lowest 
ducts occurred in warm, moist meteorological conditions (RH;85-100%, 
Ts-T:-2 to 0). 

c. Wind Speed (U) : (Figures 7, 10 and 12) 

Wind was most influential on Z* values in low wind situa- 
tions. In light winds, changes in atmospheric stability (Z/L) drasti- 
cally changed the duct. As seen in Figure 12, a change in Z/L by 
Sigma (Z/L) changed Z* from 3 to 40 meters for winds near 5 m/s. In 
moist air (RH>90%) , 1* was low regardless of wind speed. In dryer air 
(RH<85%) wind caused increased Z* as U increased from 5 to 12 m/s. 

d. Stability (Z/L): (Figures 8, 11 and 12) 

Average stability was near-neutral. The greatest deviations 
of Z* occurred when Ts-T<0 (warm air) and light wind (U<5m/s). Other- 
wise changes in Z/L, hence changes in 1* , were small as winds increased 
or colder air was advected over relatively warmer water (Ts-T>0). 

Figure 13 illustrates the interrelationships that the 
combination of selected parameters had on 1* . The ranges of Z* are 
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depicted for the primary and secondary extremes encountered in the 
joint probability density distribution. With the exception of the 
dependent parametric interrelationship of stability (Z/L) and air-sea 
tenperature (Ts-T), all parameters caused large values of 1 *. The 
primary extreme range was a reflection of the special situation in 
which large values of Z* were generated due to atmospheric conditions 
where winds were light and warm, dry, stable air overlay relatively 
cooler waters. A more significant indication of atmospheric influence 
was the secondary range. This range reflected more accurately the 
effects due to characteristic changes in atmospheric variables with 
synoptic conditions of a summer regime. 

The most influential meteorological factors were air-sea 
temperature and relative humidity, followed closely by the interrela- 
tionship of wind speed to air-sea temperature. Because the secondary 
range of wind speed to relative humidity was small, this indicated 
that large positive 1 * deviations were influenced most significantly 
by air-sea temperature. Stability only became a significant factor 
in light winds. Table IV summarizes the relationship of duct strength 
to atmospheric parameters. 
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Figure 13. Parametric Interrelationship to Duct Height 
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TABLE IV 



DUCT STRENGTH VS. AI>DSPHERIC PARANETERS 



RH C%) 



U Cm/s) Ts-T C“C) Z/L 



STRONG 

WEAK 

IVEAK 



STRONGEST 



< 80 
< 85 
85-100 
> 90 



< 5 
5-12 



< -.5 
> 0 

-2 to 0 
> 0 



> 0 



> 0 



5. Z* OCCURRENCE AND DURATION 

Histograms of duct occurrence and duration of duct heights are 
shown in Figure 14. This figiare indicates the occurrence of strong 
evaporation ducts during the summer is small and short-lived. Duct 
heights of 6m or less occurred 78% of the time. Rarely did significant 
ducts occur for more than several hours. Figure 5 shows the changing 
weather conditions throughout the summer and the effect on duct height. 
Rarely did one weather pattern dominate or persist. Consequently, the 
occurrence of strong ducting was irregular. It was most closely asso- 
ciated with atmospheric conditions that caused either cold, dry air or 
warm, dry air to move over the area. Incursions of relatively cold, 
dry air, usually from the northwest (Figinres 21- 23, 28-33 and 40-43), 
caused rapid formation of evaporation ducts (9- 12m) which rarely lasted 
longer than 9-12 hours. Dry air did not last long over the area before 
it was modified. Incursions of warm dry air were evident and produced 
brief (3-6 hrs) but strong evaporation ducts ()*8m) . Conditions less 
favorable for the formation of ducts were associated with cyclonic 
southwesterly flow of warm, moist air (Figure 20) . 
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Figure 14. Frequency of duct occurence and duration of duct 
heights . 
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Figure 15 and 16 compare the distribution frequency of 1 * . Of the 
three locations, JASIN indicates a higher percentage of low duct 
heights, Ocean Station INDIA shows its mode (group that occurs most 
frequently) in the 6-9 m range, the highest of the three. Additionally, 
INDIA had a greater percentage of Z* in the 9-18m range. Ocean Station 
.^FA located off the east coast of Greenland and in an oceanic region 
of confluence where the cold Greenland current meets the cool Irminger 
current, indicated a mode of 3- 6m and a generally higher percentage of 
1 * greater than JASI.N but less than INDIA. 

This difference between the evaporation duct climatology in the 
three areas appears to have been due to (1) the proximity of ALFA to 
the Greenland landmass and the associated, extremely cold water and 
(2) INDIA'S location in the deeper and warmer water of the North 
Atlantic current. Interpretation of visual and infrared satellite 
imagery indicates (Figures 28-33) the flow of cold, dry air off 
Greenland, coupled with cold sea surface temperatures, causes higher 
Z*' values at ALFA due to the large difference in surface humidity 
and atmospheric humidity. As the dry air mass moves across the North 
Atlantic atmospheric humidity changes are relatively small as it 
approaches both the JASIN and INDIA locations. The key difference 
between the two locations is Ts. The warmer water at INDIA causes 
higher surface humidity, thus a larger difference in atmospheric and 
surface humidity resulting in the occurrence of stronger ducts. JASIN, 
because the air mass has become more moist since leaving ALF.A, but 
enters an area of relatively waimer water, experiences the smallest 
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Figure 15. Frequency of 1 * 
occurrence (13 Jul - 6 Sep 1978) 
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Figure 16. Frequency of Z* occurrence (climatology) at 
North Atlantic Ocean Stations nearest JASIN. 
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Figure 17. Ocean Stations in the North Atlantic: ALFA 
(62N, 331V) and INDIA (58N, 19W) . 
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difference in atmospheric and surface humidities, therefore the weakest 
ducts of the three locations. 



B. HORIZONTAL NATURE OF THE EVAPORATION DUCT 

The three periods (7-8, 23-25 August and 50 August - 1 September) 
examined corresponded to times when ships were positioned at all comers 
of the JASIN triangle. The objectives of analysis for these periods 
were to: 

1. Deteimine the horizontal homogeneity of the evaporation duct. 

2. Determine which meteorological factors enhanced or degraded areal 
ducting . 

The analysis procedure for this was to: 

1. Conpare computed 2* values at different ship positions and con- 
comitant meteorological parameters to assess similarity. 

2. Analyze synoptic meteorological conditions utilizing satellite 
imagery (visual/ infrared) , ship observations and synoptic weather 
maps. 

3. Analyze the temporal variations or consistency of horizontally 
homogeneous conditions. 

4. Relate the above fields to the spatial nature of the evaporation 
duct. 

Area homogeneity was determined by comparing Z*, Ts-T, RH and U 
for ships positioned at each corner of the triangle. Azimuthal sectors 
were defined as the legs of the triangle (NW sector: NW-S, NE sector: 
NE-S, E sector: NT-NW) . Using the South (S) position as the reference, 
differences were determined relative to Northwest (NW) and Northeast 
(NE) locations (NW and NE sectors, respectively). Similarly, using the 
Northwest observer as reference, differences were determined relative 
to the Northeast location (East sector). Differences between locations 
(ships) relative to sectors are tabulated in Tables V through VII. 
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1- Z* Range in Homogeneous Atmosphere 

A necessary step to interpret the results o£ the tabulated 

differences is the performance of an error analysis to determine th'e 

possible range of computed Z* which could occur in a homogeneous 

atmospheric situation due to measurement inaccuracies. An assessment 

was required to ascertain the maximum deviation allowable in Z* for the 

three- ship system under nearly homogeneous atmospheric conditions. 

Additionally, a sensitivity analysis of I* values was performed for 

each inherent measurement uncertainty. 

Previous studies on bulk measurements [Davidson et al, 1978] 

have concluded that representative inaccuracies in bulk parameter 

measurements are ; 

RH = 5% 

Ts-T = 0.5° 

U = 10% 

Five observation times were chosen to perform the error analysis. 
These corresponded to times when air temperature, relative humidity and 
wind speed were nearly similar for all ships stationed on the triangle. 
Benchmark Z* values were computed from each of these representative 
homogeneous atmospheric situations. Atmospheric conditions at 1500 G^^^ 
August 23, 24 and 0600 GNIT 25 August were similar at all three ship 
locations. Mean values of bulk parameters were used to compute a 1* 
at each comer of the triangle assuming no measurement inaccuracies. 

A similar analysis was performed on data taken at 0600 and 1500 Qf! 31 
August. The following gives the results of the two periods: 
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CASE I 



Ts = 12 “C 


T = 


11.8 °C 


RH = 90 


SHIP 


Z(m) 


Z/L 


Z*(m) 


ENDURER 


4 


-.0055 


4.44 


MURRAY 


6 


-.0082 


4.32 


METEOR 


11 


-.0139 


4.17 



U = 8 . 1 m/ s 



CASE II 



Ts = 12.5 “C 


T = 


11.5“C 


RH = 75 


SHIP 


Z(m) 


Z/L 


Z*(m) 


ENDURER 


4 


-.015 


10.4 


MURRAY 


6 


-.023 


10.1 


^ETE0R 


11 


-.045 


9.7 



U = 10 m/s 



The 1* values computed in each of the above cases represented 
reference Z* values for each measurement height, Z (e.g. 4, 6, 11m) in 
a spatially homogeneous, mixed layer. Measurement uncertainties listed 
below were introduced into the 1* calculations and the maximum error 
computed. Percentage values for both cases reflect the change in Z* 
due only to the measurement errors of each individual parameter. 

CASE I CASE II , MEASURE^ENT ERROR (ims) 



U 

RH 

Ts 



1 % 

25% 

14% 



3 % 

8 % 



10 % 

3% 

.5 °C 



The cumulative error in the worst case situation, where measurement 
uncertainties were between two ships in a horizontally homogeneous 
atmosphere, produced the following difference in Z* values: 
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Z* (m) 



Z* Diff (Max-Min) 



4 

10 



2.7 

2.9 



Therefore, it was concluded that if the atmosphere were homo- 
geneous and the only uncertainties were in measurements of T, Ts, RH 
and U, then the greatest difference in Z* that any two ships would 
have experienced was 2.9 m. 

2. Case Studies of Horizontal Homogeneity 

Three periods were examined, in view of the error analysis, when 
ships were stationed at the comers of the JASIN triangle (approximately 
200 km apart). Simultaneous meteorological observations, taken at three- 
hour intervals, and their respective Z* values were compared. In each 
case, an average or better than average evaporation duct existed. Exten- 
sive periods of below average ducting were not considered. 

The examination included consideration of the general synoptic 
conditions, differences in local meteorological conditions between ships 
(and the factors which most influenced Z*) , the extent of the duct over 
the entire triangle (areal homogeneity) and between comers (sector homo- 
geneity). The maximum Z* difference (2.9m) between ships (due to 
measurement errors as previously described) was the criterion for duct 
homogeneity. Therefore, if the maximum difference in Z* between all 
three ships was less than or equal to 2.9m, areal homogeneity was 
assumed to exist. Similarly, if the criterion was valid between comers 
of the triangle, sector homogeneity was assumed to exist. 
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a. Case I 



(1) Period . 7 Aug (0000 GMT) - 8 Aug (2100 GJ^rr) 

(2) Observer Positions on Triangle . 

Northwest (NW) - ENDURER 
Northeast (NE) - MURRAY 
South (S) - NETEOR 

(3) Synoptic Conditions - Blocked. A ridge persisted over 
the area as it intensified and moved slowly east (Figure 19) . Northerly 
surface winds were moderate in the eastern region and lighter in the 
west. Low level instability prevailed as cold air was advected over 
warmer water. However, convection was inhibited by warm air aloft as 

an occlusion moved slowly over the area (Figures 25-27). 

(4) Results of Tabulated Differences - (Table V) . Obser- 
vations at NIV and NE were more humid than S. Winds were stronger at 
NE than NW and S. Air-sea temperattire differences were greatest at S. 
The duct heights on the average were 2m higher in the south than the 
two northern positions (i . e . 6 . 9m vs . 4 . 9m) . 

(5) Areal Homogeneity . Similar evaporation duct strengths 
occurred simultaneously at the three stations 30% of the time with the 
most sustained conditions lasting 12 hours. Figure 18 illustrates by 
sectors the percentage of homogeneous ducts and the longest duration 
experienced. Northerly flow coincided with similarity between the two 
northern positions while the greatest disparity occurred between 
southern and western positions. The most persistent horizontal homo- 
geneity occurred between NW and NE where humidities and air-sea 
teiT^erature differences were most similar. Poorest conditions existed 
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TABULAR OUFEBENCE BETWEEN SrilBS By SECTOR 
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between NW and S due primarily to the dissimilarity in Ts-T and RH and 
the consequent difference in air-sea humidity (Figures 44-47). 
b. Case II 

(1) Period : 22 August (1500 GMT) - 25 August (1500 GMT) 

(2) Observer Positions on Triangle : 

Northwest (NW) - ENDURER 
Northeast (NE) - HECLA 
South (S) - NETEOR 

(3) Synoptic Conditions - Moist Anticyclonic . An advanc- 
ing cold front passed over the area the afternoon of 23 August (Fig- 
ures 21, 28, 31). Cold, dry northwesterly air immediately followed the 
frontal passage and prevailed through 24 August (Figures 22,32, 33). 
Humidity gradually increased from the west over the area by 25 August 
as a warm front approached (Figures 23, 34, 35). 

(4) Results of Tabulated Differences -(Table VI). Generally, 
the greatest differences experienced were between N^V and NE stations 

due to the west to east movement of the air mass over the area. Winds 
and humidity were similar at all observing stations. Air-sea differ- 
ences were least at the NE position where water was the coldest, hence 
weakest ducts were generated. Duct heights ranged 1.1 -12. 2m and 
averaged 5m over the area. 

(5) Areal Homogeneity . The area experienced homogeneous 
ducting 44% of the time with the longest duration 9 hours. Homogeneity 
was nearly equal for all sectors with the northeast experiencing 
slightly stronger ducting and lasting for 24 hours following frontal 
passage (Figure 18). The greatest I* difference across the front was 
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TABULAR DUfFEREMCE BETWEEN SHIPS BY SECTOR 
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3.8m i.e. 6.6m (N^V) vs. 2.8m (NE) Horizontally homogeneous ducting 
lasted for nearly 15 hours following frontal passage before another 
impulse of dry air entered the western area. It took 15 hours for the 
dry air to move across the JASIN area and strengthen the duct in the 
east (Figures 48-55). 
c. Case III 

(1) Period . 30 August (0000 GMT) - 1 September (2100 GMT) 

(2) Observer Positions on Triangle . 

NORTHWEST (NW) - ENDURER 
NORTHEAST (NE) - HECLA 
SOUTH (S) - METEOR 

(3) Synoptic Conditions - Mobile Westerly. A weak 
depression moved over the region (Figure 19). The JASIN area was in a 
waim sector on 30 August as a cold front passed early on 31 August. A 
warm front approached from west late on 1 September (Figures 40-43.2). 

(4) Results of Tabulated Difference - (Table VII) . The 
average duct height was 6m. The most intense ducts occurred on 31 
August after the cold front passed. They occurred in the west and 
south where sea surface temperatures were highest and air-sea tempera- 
ture differences were greatest. On the average, relative humidity 
was lower to the northeast than the western and southern region which 
were largely under the influence of warm sector of the low. 

(5) Areal Homogeneity . Due to the high mobility of the 
weather systems, areal horizontal homogeneous ducting occurred only 29% 
of the time with a maximum duration of 9 hours during very low ducts. 
Sector homogeneity was poorest between the northeast and northwest 
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VARIANCE X3.5 13.0 11.6 14.2 7,8 .2 .7 .6 104.2 275.1 196 



positions and best between south and northeast (Figure 18). The maxi- 
mum difference in 1 * across the cold front was 9.5m i.e. 11.7m (MW) 
vs. 2.2m(S). The strongest duct heights associated with any weather 
system were generated on 31 August following the cold front . Though 
outside the limits of horizontal homogeneity criterion, 2.9m, duct 
heights of 8. 6-13. 4m occurred at all stations for 12 hours before 
weakening in the east (Figures 54-59). 

3. Summary of Horizontal Variation of 1 * 

The periods of the three cases studied were chosen because 
higher than normal ducts occurred for reasonably long periods . The 
horizontal nature of significant ducting in the region could be assessed. 
Each case was unique in that strong ducting was the result of differing 
synoptic meteorological conditions. 

The strongest ducts occurred as mobile westerly depressions 
moved through the area; however, horizontal homogeneity was poorest. 
Air-sea tenperature differences and relative humidity were the most 
important parameters affecting duct strength. Larger 1 * values were 
experienced at the NW and S locations because these locations had higher 
values of Ts than the .NE location. Since Ts gradients did exist between 
observation positions, sector horizontal homogeneity appeared to be 
influenced by the direction from which the air mass entered the area. 

The NW and NE positions, which experienced the largest Ts gradient, had 
more homogeneous ducting with northerly flow than westerly flow. The 
gradient between S and NE was small and resulted in consistently homo- 
geneous conditions in all three cases . 
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Case I illustrated the effect wind plays on 1 * values. Both 
iNE and-NW had high relative humidity, but NW had a greater air-sea 
temperature difference, resulting in a larger moisture flux. However 
the stronger winds in the eastern region compensated for the lower 
moisture flux producing nearly similar duct heights. A comparison of 
the S and NE stations, which had similar winds and air-sea temperature 
differences, showed significantly different duct heights. The import- 
ant parameters in this case were relative humidity and sea surface 
temperature which produced larger surface moisture fluxes. 

Case II exhibits the effects due to incursions of cold, dry air 
with moderate to fresh winds. Satellite pictures (Figures 50, 51) 
indicated low level instability under windy and drying conditions. How- 
ever, significant afternoon heating sharply reversed the air-sea tempera- 
ture difference causing the ducts to strengthen. Despite the afternoon 
variation, this case had the most consistent horizontally homogeneous 
ducting, occurring 441 of the time, as well as the most equal sector 
homogeneity of the three cases. The strongest ducts were toward the 
northwest where the air was the driest. The strong zonal flow in- 
hibited the driest from reaching the southern region where relative 
humidity varied little and air-sea temperatinre differences were small. 

Case III illustrates the effect of rapidly moving depressions 
on the formation and duration of ducts. The strongest areal ducts 
formed with the passage of a cold front but lasted only a short period 
before the next upstream weather system weakened the duct in the 
western region. The occurrence of ducts in the NE lagged the 
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occurrence of similar ducts at the NW position by six hours. 

A composite of the frequency of occurrence of areal and sector 
homogeneous ducts and maximum duration is depicted in Figure 18. It 
was concluded that areal homogeneity was a probable occurrence only 35% 
of the time for the three periods studied, with a maximum likely dura- 
tion of 12 hours. The ME sector appeared to maintain the best horizon- 
tal homogeneous conditions. The sector also appeared to exhibit good 
homogeneity but of short duration. The NW sector exhibited the least 
homogeneous conditions but when homogeneity occurred, it lasted up to 
24 hours. 
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Figure 18. Percentage of areal and sector horizontal 
homogeneity of ducting and maximum duration for three 
cases and the average. 
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VII. SUMMARY AND CONCLUSIONS 



Unique ten 5 )oral and spatial descriptions of evaporation duct 
features in the Northeast Atlantic were described from observational 
data taken from four ships spaced in a fixed array. Evaporation duct 
heights were determined from accurately measured bulk surface data 
utilizing Nfonin-Obukhov scaling in flux-profile relationships. Discern- 
ible synoptic features from surface pressure fields and DMSP satellite 
pictures were air mass histories, trajectories and horizontal variations 
(frontal zones) . Synoptic scale spatial variations of the sea surface 
teiT53erature were also available and were as important in the spatial 
variation as the atmosphere. 

Evaporation duct features established for this data set were: 

1. The mean Z* values were less than those reported in available 
climatologies based on weather ships’ observations (4.2m versus 
7.0m) . 

2. Z* values at or above 8-10 meters occur frequently in conjunction 
with outbreaks of cold, dry air from the north or northwest and 
have durations of 9 to 12 hours. 

3. Air-sea tenperature differences and relative humidity were the 
important surface layer parameters in establishing tactically 
significant values and the former was the most iirportant. 

4. If horizontal homogeneity were defined as representatively uniform 
Z* values over 200 km distances, such conditions occurred only 
35% of the times when significant duct heights existed, and for 
durations of up to 12 hours. 

Conclusions relative to synoptic scale features most closely asso- 
ciated with these near-surface features are: 

1. The primary synoptic influence on the duct height is the air-mass 
trajectory (due to its role in determining the moisture content) 
and the air- surface temperature difference. North and northwest 
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trajectories have the highest associated Z* values. 

2. The primary feature influencing horizontal homogeneity over dis- 
tances of 200km is the sea surface temperature, and its influence 
depends on the air mass trajectory. 

3. Light wind conditions modify the above specified synoptic feature- 
surface duct feature relationships. The modifications occur over 
smaller temporal and spatial scales than those of the prevailing 
features . 

A concluding assessment is that tactically significant features of 
the evaporation duct in the study area and for the study period can be 
established from synoptic scale descriptions . Further exploitation of 
available data should emphasize satellite imagery which could delineate 
surface tenperature and air mass characteristics. 

A recommendation for further efforts is that all available surface 
data from the JASIN project be used in order to extend and substantiate 
feature descriptions and relationships identified in this study. Since 
the start of this study, additional surface layer data from buoys and 
other participating ships have been compiled by the organizing institu- 
tions and made available on computer tapes . These data would increase 
the observations used by two orders of magnitude, thereby improving the 
descriptions of both the tenporal and spatial features. 
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APPENDIX A 



TABLE A- 1 

Julian (year) Day and Date 



Date 


Day 
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Day 


Date 


July 12 
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Aug. 1 


213 


Sept. 1 
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lagure 19. DAILY WEATHER MAPS 
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Figure 19. (contd) DAILY WEATHER MAPS 



00 

hs 

CK 



r> 

O 

D 

< 



>: 
— > 
00 



5 

o 

<N 




65 



Figure 19. (contd.) DAILY WEATHER MAPS 
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Figure 19. (contd.) DAILY WEATHER MAPS 
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Figure 19. (contd.) DAILY WEATHER MAPS 
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Figure 20, 500 mb (Top) and SFC Pressure Maps (Bottom) , 
1200 GMT 19 Aug 78. 
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Figure 21. 500 mb (Top) and SFC Pressure (Bottom) Maps, 1200 GMT, 



23 Aug 78. 
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Figure 22. 500 mb (Top) and SFC Pressure (Bottom) Maps, 1200 CNfl, 

24 Aug 78. 
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Figure 23. 
25 Aug 78. 



500 mb (Top) and SFC Pressure 



(Bottom) Maps, 1200 GT-rr, 
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APPENDIX B 



Defense Meteorological Satellite Program (EMSP) imagery (visual 
and infrared) of the Northeast Atlantic and high- resolution imagery 
of the JASIN area corresponding to periods of significant ducting. 




Figure 24. DNBP visual imagery, 25 July, 1212 GMT. 




Figure 25. O^ISP infrared imagery, 7 August, 1032 GMT. 
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Figure 26. EMSP visual imagery, 8 AugiiSt, 1143 GNIT. 



Figure 27. DMSP infrared imagery, 8 August, 1143 GNfT. 
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Figure 28. 



EMSP visual imagery, 22 August, 



1110 GMT. 




Figure 29. 



EMSP infrared imagery, 22 August, 1110 On. 
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Figure 30. EMSP visual imagery, 23 August, 1221 GMT. 




Figure 31. 



1>6P infrared imagery, 23 August, 1221 GNrr. 
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Figure 32. 



EMSP visual 



imagery, 24 August, 



1138 Ovfr. 




Figure 33. 



EMSP infrared imagery, 24 August, 



1138 GMT. 
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Figure 34. EMSP visual imagery, 25 August, 1054 GNfT. 




Figure 35. EMSP infrared imageiy', 25 August, 1054 ©fT. 
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Figure 36. EMSP visual imagery, 28 August, 1233 GMT. 




Figure 37. 



DMSP infrared imagery, 28 August, 1233 GMT. 
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Figure 38. 



EWSP visual imagery, 29 August, 



IISO GNfT. 




Figure 39. E^BP infrared imagery, 29 August, 1150 ©fT. 
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Figure 40. 



EMSP visual imagery, 30 August, 1107 GMT. 
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Figure 42. OVEP visual imagery, 31 August, 1218 G3yfr. 




Figure 43. I>ISP infrared imagery, 31 August, 1218 GNTT. 
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Figure 43.1 EMSP visual imagery, 1 Sept, 1135 GMT. 




Figure 43.2 



DMSP IR imagery. 



1 Sept, 1135 GMT. 
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Figure 44. High- resolution EB'ISP visual imagery and surface data. 
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Figure 45. High-resolution DNBP IR imagery 



and surface data 
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Figure 46. High- resolution DMSP \rLsual imagery and surface data 
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Figure 47. High- resolution ESVBP IR imagery and surface data. 
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Figure 48. 



High- resolution DMSP visual imagery and surface data. 
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Figure 49. High-resolution DMSP IR imagery and surface data. 
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Figure 50. High- resolution »ISP visual imagery and surface data. 
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Figure 51. High- resolution DMSP IR imagery and surface data. 
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Figure 52. 
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Figure 53. High-resolution DMSP IR imagery and surface data. 



93 





I 

I 



30 AUG (1200 mi) 



SHIP 


Z/L 


r 


IP 


T* 


a* 


u 


TS-T 


RH 


ENDURER 


>.025 


1.9 


.29 


-.04 


-.023 


7.0 


1.0 


100 


METEOR 


-.019 


2.1 


.lA 


.00 


-.017 


4.1 


.1 


95 


HECLA 


-Ml 


5.1 


.4A 


-.05 


-.055 


11.4 


1.5 


91 




Fioure 54. High- resolution DMSP visual imagery and surface data 
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Figure 5S. High-resolution DMSP IR imagery and surface data 
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31 AUG (1200GMT) 
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Figure 


56. High 


-resolution DMSP 


visual 


imagery 


and surface data. 
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Figure 57. High-resolution DMSP IR imagery and surface data. 
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Figure 58. 


High- resolution 


DMSP 


visual imagery and 


surface 


data. 
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IR imagery 


and surface data. 
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